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Abstract. We consider the interplay between excitonic and atomic motion in a 
regular, flexible chain of Rydberg atoms, extending our recent results on entanglement 
transport in Rydberg chains [Wiister et at, Phys. Rev. Lett 105 053004 (2010)]. 
In such a Rydberg chain, similar to molecular aggregates, an electronic excitation 
& ' is delocalised due to long range dipole-dipole interactions among the atoms. The 

C/3 . transport of an exciton that is initially trapped by a chain dislocation is strongly 

coupled to nuclear dynamics, forming a localised pulse of combined excitation 
and displacement. This pulse transfers entanglement between dislocated atoms 
^ (— i adiabatically along the chain. Details about the interaction and the preparation of the 

O i 1 initial state are discussed. We also present evidence that the quantum dynamics of this 

complex many-body problem can be accurately described by selected quantum-classical 
methods, which greatly simplify investigations of excitation transport in flexible chains. 

> 

' PACS numbers: 32.80.Ee, 82.20.Rp, 34.20. Cf 

(N ' 

O ■ !• Introduction 

Rydberg atoms have recently received increasing attention in cold atomic physics, to a 
large part due to their strong long-range interactions, with diverse consequences from 
dipole-blockade [H El [3] to l° n S range molecules jH |5j Ej. Among the interactions 
in cold Rydberg gases, resonant dipole-dipole interactions [7J El E] and their non- 
resonant variant (van-der-Waals interactions) [TQl [H] are particularly prominent. These 
interactions enable Rydberg ensembles to simulate the quantum dynamics of other 
long-range interacting systems, from condensed matter physics [12l [13] to molecular 
aggregates 0, [HJ [15] . We focus on the latter possibility, and explore basic consequences 
of joint dynamics of atomic motion and excitonic transport. 

Within an essential state picture, where only two Rydberg states per atom are taken 
into account, the transfer of excitation can be adequately described by using the exciton 
theory of Frenkel [16l [T7] . Following the pioneering paper by Franck and Teller [18J, 
this theory has found wide application in describing excitation transfer, e.g. in molecular 
crystals [19] , photosynthesis [20] or organic dye aggregates [21] . In all these systems the 
coupling between the exciton and nuclear degrees of freedom strongly influences the 
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excitation transfer [TH1 E21 1231 12H EH1 [23 [2H] - Similar effects will be reported in the 
present study. 

The strong interactions between the monomers of molecular aggregates lead to 
coherently delocalised entangled states [291 [301 EI] which are e.g. responsible for the J- 
band of organic dye aggregates. Recent experiments indicate robust excitonic coherence 
even in biological systems, such as photosynthetic complexes [32~ j 1331 153] . 

In all these excitonic systems the resonant nature of the interaction plays a crucial 
role. Besides the transfer of excitation, this interaction also creates a potential, which 
for an atom pair depends like 1/R 3 on their distance R. For Rydberg atoms it has been 
recognized that this potential can lead to large forces on the individual atoms [351 [T5] 
and thus cause their motion. In contrast to the atomic motion induced by van-der- 
Waals interaction, which is due to strongly off-resonant coupling, the character of the 
motion (repulsive, attractive or even mixed) in the resonant case depends strongly on 
the excitonic eigenstates [IB] . These in turn depend on the atomic positions, which 
is why excitation transport and motion become interlinked. In this respect our setup 
[T5l 136] strongly differs from that in Ref. [37], where the effect of externally enforced 
atomic motion on exciton transport is studied. 

In this article, we extend our previous studies of excitons and their dynamics 
in Rydberg chains [T5l 136] . To study exciton dynamics with Rydberg ensembles, 
one requires strong selectivity of the accessible electronic states of each atom and 
control over the initial exciton state. We consider both requirements, and furthermore 
provide additional details on the Newton's cradle type entanglement transport scenario 
reported in [36]. Specifically we vary atomic masses and interaction potentials. 
For the entanglement transport scenario, we show that two mixed quantum-classical 
methods are well suited to describe this complex many-body problem: Tullys surface 
hopping method, and the Ehrenfest method. For short chains, we validate these 
quantum-classical propagation methods by comparison with a full quantum mechanical 
calculation, finding perfect agreement. 

The paper is organized as follows: In section [2] we describe which conditions 
we imply in order to label a Rydberg chain as Rydberg aggregate. After a brief 
comparison with molecular aggregates (section F2~T1) . we describe our geometric setup and 
Hamiltonian (section 12.21) . illustrate how a simple treatment of transition dipole-dipole 
interactions can emerge (section [2731 . argue the validity of our essential states model 
(section [2 .41) . lay the basis for a description of the Rydberg chain's excitations in terms of 
excitons and their localisation (section 12. 5j) and show how the initial states for our later 
applications could be obtained (section [231) . The final part of section [2] (section [2~7l) 
details the quantum and quantum-classical formalisms used to simultaneously model 
the dynamics of atomic motion and excitons. After these preparations, we proceed 
in section H to a detailed presentation of the entanglement transport scenario first 
reported in Ref. [36] and survey the parameter space for this scenario in section HJ 
Some appendices supply further details. 
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2. Rydberg aggregates 

2.1. Brief comparison with molecular aggregates 

Since molecular aggregates have been extensively studied over the last 70 years it 
is appropriate to briefly juxtapose the Rydberg aggregates to these "conventional" 
molecular aggregates. 

Molecular aggregates appear in various contexts, ranging from organic crystals 
[T9l 138] over self-assembled cylindrical dye aggregates [39] to complex biological light 
harvesting systems [20J . These systems range from only two monomers up to thousands 
of monomers, which can aggregate into various geometrical arrangements. The 
(resonant) transition-dipole-dipole interaction between the monomers leads to entangled 
states, often accompanied by a drastic change in the absorption spectrum compared 
to that of the single monomer [211 I2H]. Besides some fundamental interest (e.g. in 
photosynthesis) the extraordinary properties of these aggregates have led to various 
applications, ranging from sensitisers in photography jlQl HI], to the measurement of 
membrane potentials [121 H3], and cancer therapy [H]. Also in the development of 
efficient, low-cost artificial light harvesting units (like organic solar cells) dye aggregates 
might play an important role [4"5] Ef6] . 

In molecular aggregates the monomers are held at their positions and orientations 
e.g. by a protein environment or by van-der-Waals interactions, with distances of the 
order of a few Angstrom. In the Rydberg aggregates investigated here, the distances 
are of the order of a few micrometers and, most importantly, the Rydberg atoms are 
free to move. The main difference, however, between molecular aggregates and Rydberg 
aggregates is the internal structure and the environment. While the Rydberg atoms are 
at ultra-cold temperatures and interact only weakly with the environment, the electronic 
excitation in the molecular case does strongly couple to the environment (often at 
ambient temperature) and a plethora of internal vibrational modes [4"7j . This typically 
necessitates various approximations and assumptions in the theoretical description of 
molecular aggregates , since often details about the environment or even the precise 
arrangement of the monomers are unknown. Furthermore, due to the small distances, 
the direct experimental observation of coherent energy transfer in molecular aggregates 
is challenging. Hence, related investigations are typically of spectroscopic nature [48J, 
and infer the exciton dynamics only indirectly. 

In contrast the beauty of Rydberg aggregates is, that individual excitation and 
manipulation of the atoms can be done more easily. Also, since environment and 
vibrations do not play a role, it is possible to develop a detailed theoretical description 
where common approximations can be checked. 

2.2. General setup 

We study a chain of N identical atoms with mass M and denote the position of 
the nth atom by R n . These positions are grouped into a 3N dimensional vector 



Adiabatic entanglement transport in Rydberg aggregates 



4 



|3 5p> 

|35s> 




I 2s > _ 

atom 1 2 3 4 5 



Figure 1. Visualization of the electronic state | 7Tl ). 

R = (R 1; . . . , R^v) 7 !]]. In the following we will refer to these coordinates as nuclear 
coordinates. Each atom should be initially well localised, for example in the ground 
state of an optical lattice or a micro- lens array [49]. We can then ensure that the 
distance R nm = |R n — R m | between the atoms is large enough to neglect the overlap 
between the electronic wave functions of atoms n and m. 

Consider a situation where all but one of the N atoms are in a Rydberg state 
with principal quantum number v and angular momentum 1 = 0. The remaining atom 
is in an angular momentum I = 1 state \vp), which we will call the "excited" state. We 
now define the single-excitation Hilbert space, whose electronic part is spanned by 

|7T re ) = \s---p---s), (1) 

a state in which atom n is in the p state and all others are in the s state. In figure [1] 
the state 1 7Ti ) and the initial spatial arrangement is sketched for the case N = 5. Note 
that the role of s and p is more or less interchangeable. 

For a clear cut picture of exciton transport, the states (JT|) should form the essential 
part of the electronic basis for the whole aggregate. This requires that transitions to 
other states, such as \ s • • -p - • • s • • -p - • • s) or | s • ■ ■ d • ■ ■ s ) are negligible, because they 
are energetically far detuned with respect to the relevant couplings. Then, the only 
relevant interactions occur within the space spanned by ([I]) and conserve the number 
of excitations. In section 12.41 we show by an example that these requirements can be 
fulfilled. 

In terms of the basis (JT|), our total Hamiltonian describing atomic motion and 
interactions within the essential states manifold is given in atomic units by 

ff (R) = -Ew +Fl(R) ' (2) 

n=l 

Here, the electronic Hamiltonian, which depends on the nuclear coordinates, is 

# Cl (R) = Vnm(Rnm) | 7T re > ( 7T m | , (3) 



| We present our introductory theory as far as possible in three dimensions (3D). For all our results we 
only consider one-dimensional motion (ID), assuming transverse motion is frozen out by the confinement 
of the atoms. 
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where 

V nm (Rnm) = (-ir-lf- ( 4 ) 

is the dipole-dipole coupling between atoms n and m and -R nm their separation. We 
parameterised the strength of the coupling by its magnitude /i 2 and sign 77 G {0,1}. 
Due to this resonant dipole-dipole interaction the "excitation" can be transferred from 
atom n to m. We outline in section 12.31 why we can avoid a more complicated, angular 
dependent [9] expression. 

For most specific examples throughout this article, we consider the atomic species 
7 Li Lithium. Among the work horses of cold atom physics, this atom is one of the 
lightest and hence most suited to display phenomena of dipole-dipole interaction induced 
motion, within the time-scales available. Its atomic mass is roughly M = 11000 au and 
transition dipole moment have the strength /x = 1000 au between s and p states with a 
principal quantum number v « 30 . . . 40. 

In section H] we survey the response of dynamics dictated by ([2]) to changes 
of Hamiltonian parameters. To this end we also generalise the type of interaction, 
considering V nm (R nm ) = (— l^/i 2 //?^, where a can for example vary from a = 1 ... 6, 
with character of the interaction potentials ranging from Coulombic to van-der-Waals. 
Keep in mind though, that unlike conventional Coulomb or van-der-Waals interactions, 
those considered here would still have a resonant transition character. 



2.3. Dipole-dipole interactions 

In this section we outline how the simple form V nm (R nm ) = (— l) n ii 2 /R^ im can be 
obtained for dipole-dipole interactions. For this purpose, we consider a binary atom 
system with separation R nm = R m — R n and define R = |R nm | and R = H nm /R. We 
assume one of the atoms is in a \vs) state and the other in a \vp) state, where v is 
the (large) principle quantum number, subsequently suppressed. As long as one ignores 
directional effects, the essential states Hilbert space for such two atoms is spanned by 
\sp), \ps). Considering angular dependent transition dipole-dipole interactions amounts 
to taking into account also the magnetic quantum number. We then have six essential 
states: \{p,l}s), \{p, 0}s), \{p,—l}s), \s{p, 1}), \s{p, 0}), \s{p, — 1}), using an obvious 
notation that writes the magnetic quantum number m of the atom with / = 1 within 
the curly brackets. 

The non-vanishing dipole-dipole transition amplitudes between those states are j9] 

Vlm,00]00,lm' ~ ~ \j ~ K 1 ) 




2 



-m' m! — m 



*V'-m(R.), (5) 



where Yj jm are spherical harmonics and (• • •) denotes the Wigner 3j coefficient. The 
matrix-element Vi imi) i 2m2 .i> im > it i> 2m > 2 describes a transition between the two-atom states 
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indicated with primed and non-primed subscripts. d v \,vQ is the radial overlap-matrix 
element between the I = and I = 1 states. We refer to Ref. j9] for further details. In 
matrix-form, using the above basis ordering, one obtains 



V 



with sub-matrices 



V, 

' ps 



ps 





(6) 



ps 



R? 



Scos^ 



6 



4= cos 9 sin 9 



^y=- cos 9 sin 9 

1-3 cos 2 e 



cos 9 sin 9 



■ cos 9 sin 9 



V2 
3 cos 2 0-1 



(7) 



In this matrix, the element (V ps )ij contains the amplitude of transitions from a state 
\{p,rrij}s) to \s{p, mi})), where mi,mj e {1,0,-1}. We used the short-hand jl 2 = 
{dn a i,n b o) 2 ■ The angles 9 and <p describe R in a spherical polar co-ordinate system 
whose z-axis (z) is the quantization axis with respect to which the magnetic quantum 
numbers m is defined. A useful choice of z will be given by the polarisation direction of 
the light-field used for the initial-state creation, see section 12.61 

We will consider two specific simple cases, assuming a linear Rydberg chain. 
case (i): Choose z along the direction of the chain. Then for all distance vectors R nm 
we have 9 = and 



V ps 



3i? 3 



1 
0-2 
1 



(8) 



Thus the magnetic quantum number of the excitation is conserved. Depending on the 
selected magnetic quantum number m, we can realise different signs 77 and magnitudes 
of the interaction. 

case (ii): Choose z perpendicular to the direction of the chain, which we assume to be 
in the x direction. We then have (setting 9 = 71/ 2 and = 0) 



ps 



6i? 3 





3 




2 




3 




1 



(9) 



It can be seen that the m = state decouples and yields a dipole-dipole interaction 
transport without angular-dependence. For all choices of quantisation axis and magnetic 
quantum state, we finally define the parameter fi 2 used in section 12.21 as the modulus of 
the factor multiplying R~ 3 . How a specific magnetic quantum- number for the excitation 
can be realised is described in section 12.61 



2.4- Validity of essential states model 

Dipole-dipole transitions in principle do not only take place within the | sp ) , | ps ) 
manifold, but also to other states not included in our essential states picture [9]. 
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Nonetheless, parameters where these other states can be ignored can easily be found, 
as we demonstrate now. 

We exploit that transition probabilities from the two original states to all other 
di-atomic states are negligibly small due to their energy mismatch or selection rules. In 
other words, all two-atom states to which a direct dipole-dipole transition is possible 
are much farther detuned than the strength of the transition matrix element. We will 
illustrate this argument in the following for 7 Li. 

Consider the two-atom states |35,s) \35,p), | 35, p ) ® |35,s). The two-atom 
states energetically nearest and connected via a single dipole-dipole transition are 
| 36, s) <S> | 34, p), | 34, p) ® | 36, s), detuned from our essential states by A = 8.78 GHz 
and connected with a coupling strength of about V = 65 MHz^jl A simple analytical 
four state model for detuned Rabi-oscillations then predicts population transfer out of 
the essential-state system of the order of (V/A) 2 = 5 x 10 -5 for time scales considered 
in this paper. 

As more rigorous justification of the essential state mode, also accounting for 
successive, cascaded transitions, we propagated the state | 35, s) | 35, p) + | 35, p ) ® 
| 35, s ) (an exciton eigenstate, see section 1231) with a Hamiltonian that contains all u,l,m 
states from v = 34 to v = 36, setting all dipole-allowed transition matrix elements V 
to V = fi 2 /d 3 for fi = 1000 au and d = 2/im. This value for /z overestimates almost all 
transition dipoles and the value for d is the smallest separation occurring in the atomic 
dynamics of section [3) Within this "worst-case" scenario, total transitions out of our 
target essential states manifold are of the order of 1 x 10~ 4 within 20/is, which is longer 
than the simulated time-span in section [31 

Finally kinetic and potential energies of the dynamics presented in our manuscript 
amount to only about 5% of the energetic separation between | 35p), | 35d). 

These estimates, while exemplary, show that there is no general problem with 
finding physically realistic scenarios with Rydberg atoms that can be described well 
with our model. 

2. 5. Excitons, exciton localisation and full aggregate initial state 

To gain some insight into the structure of the dynamics induced by the Hamiltonian 
(T5]), consider eigenstates of the electronic Hamiltonian 



For each R there are N eigenstates labeled by the index k. Each of these eigenstates 
can be expanded in terms of the previously introduced basis | vr n ) as 



These eigenstates are termed Frenkel "excitons" [HI [T7j and form an adiabatic 
(Born-Oppenheimer) basis in the language of molecular physics. The corresponding 

§ Wc obtain the lithium state energies using quantum-defect theory as described in |50j and transition 
strength as outlined in [9] with numerical Numerov calculation of radial overlaps. 



iT(R)|v? fc (R)> = f4(R)|<MR)>. 



(10) 




(11) 



■111 
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Figure 2. Initial population of the (out of phase) repulsive excitonic state for different 
a as a function of the dislocation ratio a/xo- The vertical black line indicates our 
choices for a/xo employed in section [4] Other lines indicate the population in the 
different excitonic states for a chain containing 6 Rydberg atoms. 




Xq 



Figure 3. Sketch of the initial total density distribution of Rydberg atoms for the 
case of N = 5 atoms. 



eigenenergies [/^(R), which also depend parametrically on the nuclear coordinates R, 
define the adiabatic potential surfaces. As evident from (jlip . an exciton is a coherent 
superposition of different localised excitation states. 

Consider now a regular chain of Rydberg atoms with spacing xq, which is perturbed 
by a dislocation of two atoms in close mutual proximity (distance a), see figure The 
interaction between these atoms is much larger than interactions involving the remaining 
atoms. As a consequence two of the exciton states localise on the dislocation atoms. 
For a < io the state whose Born-Opp. surface has repulsive character [15] can be 
approximately written as | </? rep ) « (| 7Ti ) + (-1)^1 tt 2 ))/ v / 2. Such repul sive dimer states 
are observed e.g. in [8]. In figure |2] the excitonic population on the various atoms as 
a function of a/x^ is shown for the case N = 6 and different interaction exponents 
a. Our survey of dynamics presented later requires sufficiently good coherent exciton 
localisation on the dislocation atoms, which lead to our choices of a/xo indicated by the 
vertical black lines in the figure. 
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We now are in position to discuss the whole initial state of our Rydberg aggregate, 
describing the electronic state and the position of the atoms. The initial spatial wave 
function of each atom is assumed Gaussian with standard deviation ctq. This resembles 
an experimental situation where the ground-state atoms are trapped in harmonic 
potentials prior to their excitation to the Rydberg level, as discussed in the next section. 
For a sketch of this arrangement, see figure [3l 

We take the complete initial wave function (i.e. containing nuclear and excitonic 
degrees of freedom) as 

TV 

| *(t = 0) ) = | <^ rep (R) ) J] G (R n ), (12) 

n=l 

G (R n ) =Afexp(-|R„-R 0ri |72<To)> (13) 

where Ro n is the center of mass of the n-th Gaussian and M a normalization factor. 
As above, the index "rep" denotes the unique exciton state with globally repulsive 
behaviour [15]. In (fT2"|) . the exciton state varies with the nuclear co-ordinates R within 
the initial Gaussian distribution, but as long as do C a this effect will be small. We 
found that one obtains almost identical dynamics to that presented in the following, if 
<^ rep (R) is replaced by (^ rep (Ro) in (TT2]) . 

2.6. State preparation 

In this section we will briefly describe how states discussed in the previous section could 
be prepared experimentally. 

At the very beginning, ground state atoms are confined e.g. in optical traps created 
by micro lens arrays that provide the desired spacial arrangement, in our case a linear 
chain with distances a and xq. Ideally there would be exactly one atom per trap-site. 
These ground state atoms are then transferred into a certain Rydberg state (say 35s) 
via laser excitatiorjjj]. We now have a state, which we denote by |s---s), where all 
atoms are in the same Rydberg state | 35s ). Due to the ultra-cold temperatures the 
Rydberg atoms can be regarded as frozen (the distances between the sites have to be 
chosen such that acceleration and blockade effects due to van-der-Waals interaction 
[5T] are negligible). Then by applying a short microwave pulse, which will be specified 
below, the eigenstate of the chain | ^ ini } « (| ni ) + \ ix 2 ))/v / 2 can be excited. To reach 
I ^mi ) ~ (| tti ) — 1 7T2 }) / ' V% we require a further phase-flip described below. 

For definiteness we choose the propagation direction of the microwave pulse 
perpendicular to the chain. Within the dipole approximation the interaction of atom n 
with the microwave pulse is given by 

W n (t) = -» n E(t) (14) 

with the dipole operator \i n and electric field £(t). 

|| If there is more than one atom in the ground state trap one can ensure, due to dipole-blockade, that 
there is actually only one Rydberg excitation per site. 
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Since our target initial electronic state is essentially located on two atoms, in the 
following we discuss the microwave excitation exemplarily for a dimer. The extension to 
larger systems can be easily done. In this section, we enlarge our essential state space 
beyond 1 7Ti ) and 1 7r 2 ) to include also the "ground state" | ss ) and the doubly excited 
state | {p, mi}{p, 777-2} )• In this basis the Hamiltonian of the system can be written as 



H 



2E S 


W 1 (t) 


W 3 (t) 





W\(t) 


(E s + E p )l 3 




W 2 (t) 


wl(t) 


v ps 


(E s + E p )l 3 


w^t) 







W\{t) 


2E, p lg 



(15) 



with Vp S given by ([7]) and E s and E p denoting the energies of the respective Rydberg 
states. Further, l n denotes a n x n unit-matrix, Wj(t) a 3 x 1 vector and Wj (t) a 3 x 9 
matrix. The components m of Wj(t) are given by (Wj(t)) m = (ss \Wj(t)\ {p,m}s), and 
similarly W_j(t) has elements given by ( s{p, m 2 } \W\(t)\ {p, mi}{p, m 2 } ) 

We now take the microwave to be linearly polarized and choose our quantization 
axis z in the direction of the polarization, i.e. £(t) = S(t)z. Similar as in Section 12731 we 
get for the matrix elements 

d U Q,ul r 



( ss \ fJ>iz\ {p, m}s ) 



JmO 



(16) 



From the results of section I2.3[ we then notice that microwave polarisation along the 
chain (see (IE])) or perpendicular to the chain (see (|9])) leads to a de-coupling of the 
angular momentum state m = from the other m states. Thus we can consider the 
reduced Hamiltonian 



H 



2E S 




n 2 (t) 







{E a + E p ) 


Vps 


n 2 (*) 


n 2 (t) 


V ps 


(E s + E p ) 







n 2 {t) 




2E P 



(17) 



with Q n (t) = -j$E{t)d ^0,1/1 and VL n is of the order of Vt n . It is instructive to diagonalize 
within the one-exciton space to obtain the "eigenstates" | ± ) = -^(| i\\ ) ± | tt 2 )) with 
energies E± = (E, 
in this basis 



Ep)±V ps . Taking E s as zero of energy we obtain for the Hamiltonian 



H 





n + (t) 
n_(t) 




n + (t) 

E p + Vps 


n + (t) 





E p -V % 



ps 





n + (t) 
n_(t) 

2 En 



with n ± = -^(Qi ± n 2 ). Since we are dealing with identical atoms = 0, hence the 
microwave couples only to the symmetric state | + ) = -^(l^i) + I n 2 ))• Thus in order 
to be resonant with the transition \ ss) — > | + ) we will detune the microwave by V ps 
w.r.t. the atomic transition frequency. This also means that the microwave is detuned 
by 2V ps w.r.t the transition from the state | + ) to the doubly excited state | pp ), so that 
the population of the doubly excited state will be strongly suppressed. 



Adiabatic entanglement transport in Rydberg aggregates 



11 



Ideally the microwave pulse should transfer all the population from the | ss ) state 
to the | + ) state and be so short that the atoms do not move appreciably during the 
duration of the pulse. We have done full numerical simulations of this excitation scheme 
for three Lithium atoms, and we found that pulses of few nanoseconds duration can be 
used to achieve this goal. 

It also is of interest to access the aggregate eigenstate | — ) . We will actually 
focus on dynamics arising from a | — )-type initialstate and interactions with 77 = 1 
throughout this article, since in that case the smaller energetic separation between the 
totally repulsive adiabatic state and its energetic neighbour leads to more interesting 
non-adiabatic effects. This scenario was also considered in our previous work [T5], [36] • 

Since | — } does not directly couple to the linear polarised microwave, as argued 
above, this requires a second state preparation step in which e.g. the phase of the 1 7r 2 ) 
component of the quantum state is inverted. This can be achieved using a Rabi-27r laser 
pulse, which is resonant on the transition from | up) to e.g. the absolute ground-state 
I 2s) and spatially focussed to only interact with atom two [52] . 

2. 7. Dynamical methods 

Up to this point we have introduced the Rydberg aggregate as an ensemble of alkali 
atoms with parameters chosen to enable a description of collective excitations in terms 
of Frenkel excitons, and explained how the atoms can be brought into the required 
internal electronic states. To form a flexible Rydberg aggregate, we further wish to 
include motion of the atoms. We now list different possibilities to describe this motion 
numerically. 

2.7.1. Exact solution: Schrddinger's equation The full quantum-mechanical many- 
body problem posed by the Hamiltonian ([2]) is conceptually straightforward, but 
becomes quickly intractable as the number of atoms N is increased. For small N it 
is however no problem to directly solve the Schrodinger equation 



Expanding the full wave function in electronic (diabatic) states according to | ^(R) } = 
Hn=i 0n(R)| vr n ), we arrive at 



We solve (1201) for three Li atoms in order to validate the quantum-classical methods 
presented further below, which in turn will then be faithfully used for longer chains. In 
practise, the irrelevant centre-of-mass degree of freedom is removed from (120]) resulting 
in an effectively two-dimensional (2D) problem. This is solved on a discrete spatial grid. 
The above diabatic representation of the wave function </> n (R) allows a straight 



(19) 




(20) 
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forward propagation. To interpret the results and compare them with the quantum- 
classical methods, it can also be beneficial to move to the adiabatic representation 

N 

|v&(R)) = ^0 fc (R)|^(R)). (21) 

k=l 

The two representations are related by 

fe (R) = J]O fen (R)0 n (R). (22) 

n 

with O kn (R) = ( <£>fc(R) 1 7r re ). For instance, the initial state (fT2l) corresponds to 
rcp (R) = n^Li 0c(-Rn) and 0fc(R) = for k ^ rep in this representation. 

When analysing our results, we will not show the full N- dimensional nuclear/ atomic 
wave function but focus on the more intuitive total atomic density, which is given by 

n n r 

n ^ = E E / d^R^mm 2 . (23) 

j=l m=l 

Here J d N ~ 1 H{jy denotes integration over all but the jth nuclear coordinate. The density 
n(R) gives the probability to find an atom at position R. 

We will assume that wave functions of different atoms never occupy the same space. 
For the calculations shown, this assumption turned out to be valid. 

2.7.2. Quantum- classical propagation When the number of atoms N exceeds values 
where the direct quantum solution of the time-dependent Schrodinger equation (1201) is 
tractable, we resort to mixed quantum-classical methods, namely the Ehrenfest method 
[531 l5~4] and Tully's fewest switching algorithm jSU [55]. In both approaches, the 
nuclear coordinates R are treated classically and an ensemble of trajectories R(£) is 
propagated in a way specified below. In order to represent the initial nuclear wave 
packet, we randomise the initial positions and velocities for the trajectories according to 
the Wigner distribution of the initial state (1121) . Since the spatial density of each atom 
is assumed to be Gaussian, this simply amounts to un-correlated Gaussian spread of 
both, position (with standard deviation cr /v^2) and velocities (with standard deviation 
l/(\/2croM)). To obtain the total atomic density n{R), the positions of the atoms are 
binned throughout all trajectories. 

The excitonic propagation is done by expanding | ^(R, t) ) = J2k=i^k(t)\ ^(R) ), 
where the complex amplitudes c k are determined by 

d N 
i— c k = U k (R)c k - i^2~R ■ d kq c g , (24) 

q=l 

where £4(R) are the adiabatic potential energy surfaces defined in (TTTJT) and 

d kq = (^(R)|V R |^(R)) (25) 
are the so-called non-adiabatic coupling vectors. 
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The two methods differ in the classical propagation method for the nuclear 
coordinates. In the Ehrenfest method the nuclear dynamics is determined by Newton's 
equations 

MR = - V R C7(R, t) (26) 

with the average potential U(R,t) = ( *(R, t) |# cl (R)| *(R, t) ) = J2 k (*) \ 2 U k (R). 

In contrast in Tully's method each trajectory moves classically on a single adiabatic 
surface £4(R), except for the possibility of instantaneous jumps among the adiabatic 
states. Between jumps the classical equation of motion is 

MR=-V R f4(R). (27) 

Details on Tully's method and our numerical implementation are given in Appendix B| 
Now, when performing the average over trajectories the spreading due to the surface 
hopping is combined with the spreading due to different trajectories for different initial 
classical nuclear positions. 



3. Entanglement transport 

In the previous section we have explained the design of a flexible Rydberg aggregate and 
our various methods for dynamical propagation. Hence we are ready to consider the 
dynamical problem introduced in Rcf. [36J in more detail. We study the effect of resonant 
dipole-dipole interactions on a regular linear chain of Rydberg atoms. Initially we impose 
a "deformation" in the distances between the atoms that gives rise to an associated 
localised exciton state, which is strongly repulsive. We demonstrate a strong correlation 
between the resulting exciton dynamics and the motion of the atoms. A combined pulse 
of atomic displacements ("deformation") and localised electronic excitation propagates 
adiabatically through the chain in a manner reminiscent of Newton's cradle. We show 
that this can also be viewed as adiabatic entanglement transport, since the initial 
electronic state | <^ rep (R, t = 0) ) is a Bell-state [56]. To see this, we re-write the initial 
state 

|Vxep(R,* = 0)>«-^(|7r 1 )-|7r 2 >) (28) 

<8>|s---s), (29) 

where the state in square brackets concerns the dislocated atoms, and | s ■ • • s ) the rest 
of the chain. Prior to demonstrating the combined transport of displacement, excitation 
and entanglement, we validate the quantum-classical methods required for larger chains. 



i 

7/5 



ps) — I sp) 



3.1. Comparison and validation of methods 

To confirm the applicability of quantum-classical numerical treatments to the dislocated 
chain of section 12. 2[ we consider the smallest nontrivial chain, namely A^ = 3. In this 
case it is no problem to solve the full Schrodinger equation numerically exactly. We are 
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then in a position to compare all three propagation schemes outlined in section 12 .7\ full 
quantum-mechanics (QM), Tully's fewest switching (Tully) and the Ehrenfest method 
(EF). We consider two distinct scenarios: (i) predominantly adiabatic dynamics for the 
validation of the quantum-classical methods for the subsequent section 13.21 (ii) strongly 
non-adiabatic dynamics, in order to highlight the differences in propagation algorithms. 

Scenario (i) is shown in figure @Ja-e). We used M = 11000 au and \i = 1000 au as 
in Ref. [36]. The quantum mechanical probability to find an atom at a certain position 
predicted by QM and the corresponding semiclassical methods show perfect agreement. 
As the dynamics is almost completely adiabatic, each avoided collision between two 
atoms is accompanied by excitation transfer. We will highlight this in detail in the 
next section, where we consider longer chains. Note that Tully and EF even perfectly 
reproduce the small fraction of population that has switched to the neighbouring surface 
0, as can be seen in Fig. @] (b). 

For scenario (ii), shown in figure H^f-j ) , we changed our parameters to M — 1800 
(hydrogen) and \x = 200 au. This fictitious scenario was solely chosen to increase the 
system's non-adiabaticity and is probably not realistic. Due to increased diffusion and 
collisions, we extended all three models by a phenomenological treatment of ionisation, 
presented and justified in Appendix C It can be seen that in contrast to scenario (i), 
there are now significant transitions from the initial surface "rep" to "mid". On this 
surface, the trimer no longer feels an overall repulsive potential [15]. Consequently, 
atoms that have undergone a change of adiabatic state can approach each other closely 
where they ionise. This is reflected in the drop of overall population for the QM and 
Tully models. In contrast, atoms in the EF model always propagate according to a 
state averaged potential, which due to 75 % population on the repulsive surface is still 
dominantly repulsive. Consequently we do not observe significant ionisation in the EF 
model. Despite this main difference, it can be seen that the overall state-population as 
well as spatial density distribution of the exact QM model is still fairly well reproduced 
by both Tully and EF. 

The physical situation shown in figure HJ^a-e) is quite similar to that presented in 
the following section, except for the number of atoms. The quality of agreement between 
the three disparate methods found in the case N = 3 gives confidence that the quantum 
classical methods will produce reliable results also for the longer chain considered next, 
for which a solution of the Schrodinger equation would no longer be feasible. 



3.2. Coupled atomic and electronic dynamics 

The atomic motion and excitation transfer for a chain of iV = 7 atoms, when starting in 
the exciton state with highest energy (which corresponds to the fully repulsive state) is 
shown in figure |5] and figure EJ As expected, initially the two close atoms strongly repel 

% For the case N = 3 there are three adiabatic surfaces, one overall repulsive ("rep"), one attractive 
( "att" ) and one energetically between those, which we label "mid" . 
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Figure 4. Comparison of QM, Tully and EF for nuclear dynamics in the case 
N = 3 for two different parameter sets, yielding adiabatic (a-c) and non-adiabatic 
dynamics (f-j). The time evolution of the total atomic density n(R, t) (a,f) is shown 
together with a comparison of Tully's surface hopping calculations (black solid line) 
with the full quantum evolution (red dashed line) and EF (blue dotted line) in the 
other panels, (b) Relative population on the energetically nearest adiabatic surface, 
n mid = / ^R-l0mid(R)| 2 ( n mid = l^micd 2 m the Tully /EF algorithms), as a measure 
of the propensity of non-adiabatic transitions. The index "mid" is defined in the text. 
The inset shows the differences Tully-QM (black solid line) and EF-QM (red dashed 
line), (g) Similar to (b) but showing the total population (Tully (blue solid line), QM 
(green dashed line), EF (blue dotted line)) and population of all three surfaces "rep", 
"mid", "att" in descending order. (c,h) Spatial slice n(x, ti), with t\ as indicated by 
the first vertical white lines in (a,f). (d,i) Spatial slice n(x, t^), with ti as indicated by 
the second vertical white lines in (a,f). (e,j) Difference Tully-QM and EF-QM for the 
density profiles at t 2 with lines as in the inset of (b). 



each other. When atom 2 approaches atom 3 these atoms start to repel each other. 
Atom 2 slows down and atom 3 accelerates. In this way the momentum of atom 2 is 
transferred through the chain to atom 7, which is reached at t « 5.5/is. Then atom 7 
moves away from the remaining N — 2 atoms, as atom 1 did already at the beginning of 
the evolution. The remaining atoms form a regular chain with distance xq between the 
atoms and positions shifted by xq — a w.r.t. the initial position of the respective atom. 
This chain is in a repulsive state and the atoms drift very slowly apart in a manner 
typical for a regular chain as discussed in [T5] . 

Note the spreading of the initially quite localised wavepackets right from the 
start, for example atom 1. This is due to the initial spatial width of the Rydberg 
atom distribution cx , which gets converted into strong velocity spread Av due to the 
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Figure 5. Dynamics of atomic motion and excitation transfer. Shown are snap-shots 
in time as labelled, of the total atomic density, as a function of spatial co-ordinate. 
The colour shading reflects the probability of the underlying atom to be excited, and 
hence demonstrates exciton transport. 



steep slope of the dipole-dipole potential. One expects Av = (2<7o/U 2 /Ma 4 ) 1//2 . Then 
estimating Ax = Avt describes the spreading of atom 1 well. Atom 2 initially obtains the 
same large velocity spread, in the following elastic collision this is however exchanged 
completely against the (narrow) velocity distribution of 3. After the dislocation has 
traversed the chain, only the outer atoms have a considerable spread in velocity which 
results in a large position smear as time progresses (see figure [5]). During this transfer 
of momentum there is negligible overlap of the spatial distributions of different atoms, 
even at the avoided collisionaE). Up to now we have restricted our discussion to 
the "slow" , "macroscopic" movement of the Rydberg atoms. The interaction strength 
between a pair of Rydberg atoms at a distance a = 3.5/iin (corresponding to the average 
closest distance between atoms up to the time t=5.5//s) is approximately 141MHz, which 
corresponds to a "transfer" time of 0.02 /xs. This is much faster than the time scale 
of motion of the atoms. The colour shading in figure [3J and figure [6b visualise how 
the electronic excitation evolves in time. One sees that the excitation gets transferred 
such that it is always localised on the two instantaneously nearest atoms, in accordance 
with the structure of exciton eigenstates outlined in [13] . After 5.5 /is the momentum 
that was transferred through the chain kicks out the last atom, and a well defined 

+ In interpreting figure [5] keep in mind that a narrow gap between the total density peaks associated 
with two neighbouring atoms does not imply that the atoms approach closely: correlations between 
atomic positions are strong and result in the absence of actual close encounters. 
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Figure 6. Dynamics of atomic motion and excitation transfer averaged over 10 5 
realisations, (a) Total atomic density as denned in (|^|) . obtained by binning the 
classical trajectories R(t), see | Appendix B| . (b) Diabatic populations |c m (£)| 2 in 
the localised state basis \^(t)} = J2 n c n (t)\ ir n ). The row m shows the excitation 
probability of atom m . (c) Population on the adiabatic (eigen-) surface rep (red) and 
the energetically next one (blue), (d) Binary entanglement (see |Appcndix A| E n:Tl+ i 
for neighbouring atoms. The pair of indices n, n + 1 is assigned to each line near its 
maximum. 



close proximity pair no longer exists. The exciton state then assumes the shape for 
an equidistant chain, de-localised over the entire chain (consisting of the remaining N- 
2 atoms), which subsequently slowly spreads out. However this state change is not 
completely adiabatic as can be seen in figure Eb where the adiabatic population on the 
initial (repulsive) adiabatic surface together with the population on the neighbouring 
adiabatic surface is shown. One clearly sees a change of population around t = 5.5/is, 
which is the time when atom 7 starts to separate from the chain. The duration over 
which population transfer between the surfaces occurs corresponds to the time during 
which the excitation localised on atom 6 and 7 spreads over the remaining chain (see 
figure [5] and figure |6b). This change of the adiabatic populations can be understood 
in a simple way: As noted above, up to t « 5/is the excitonic transfer time was much 
faster than the nuclear dynamics. During the time in which atom 7 leaves the chain, 
however, the excitation has to delocalise over the whole remaining chain to stay in the 
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fully repulsive adiabatic state. The distances involved in this redistribution of excitation 
are much larger than a or i , hence the electronic time scale is slower, now of the order 
of the nuclear motion of atom 7. After the delocalised state is reached, the relevant 
nuclear dynamics becomes very slow - the system behaves adiabatic again. We found 
that the magnitude of these non-adiabatic transitions increases with chain length N if 
all other parameters are kept constant, reflecting a decrease in the energetic separation 
of the involved adiabatic states for larger N. 

So far we have viewed the dynamics of excitation transport essentially as a wave 
spreading phenomenon on a chain whose constituents are free to move. It is possible 
to give the observed phenomenon a quite different twist, by considering the dynamical 
transport of entanglement that is linked to the excitation migration. In particular we 
focus on entanglement within the subsystem comprised of the electronic state of atoms 
n and m only, denoted by E n>m . This subsystem can contain much less information 
than the full many-body quantum state, hence entanglement therein is expected to be 



more robust and simultaneously more accessible. We summarize in |Appendix A| how 
we calculate the relevant bipartite entanglement of formation [57]. As can be seen in 
Fig. [6] (d), the initially perfect entanglement between 1 and 2 is transported through 
the chain with only minor losses up to the point where the final atom leaves the chain Q 
Then the exciton state de-localizes over the entire chain, with a resulting drop of bipartite 
entanglement. A comparison of panels (b) and (d) of figure [6] makes it apparent that 
entanglement is here a direct consequence of coherent, delocalised excitation: Whenever 
the diabatic population on both members of a neigbouring atom pair is large, so is the 
mutual entanglement. 

4. Parameter dependence of the entanglement transport 

In the following we will investigate how the coupled excitonic and nuclear dynamics 
depends on the mass M of the atom, the magnitude [i of the used transition dipole and 
on the absolute initial positions of the atoms. These dependencies can be in principle 
be studied in experiment. In addition we will also investigate changes of the functional 
form of the long-range interaction, which is more of general theoretical interest. As in 
the previous section, for the following calculations we will use and compare the two 
mixed quantum classical methods. 

Since we focus on dynamics which is more or less adiabatic, the motion of the atoms 
is approximately governed by 

£)2 2 2 

M—R = -V R U rep (R) = -V R ^ = 3^, (30) 

which is equivalent to J^-R = ^jj^i ■ If we now scale 

Rx = XR, (31) 

* We have verified that both methods outlined previously give the same entanglement evolution for 
the case N = 3. 
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Figure 7. Entanglement E nm (t) between neighbouring atoms n and m for different 
values of the transition dipole moment \x calculated with Tully and Ehrenfest. For each 
[i the time is scaled with to = io(^) = T/n with T = 6.44/is. The curves for different 
fi are indistiguishable. 



M p = pM, (32) 

lh =7^» ( 33 ) 
we see that (130|) remains invariant, if time is also scaled by 

r = ty-^W 2 . (34) 

This means, for example, that for doubled transition dipole moment \x one expects the 
dynamics to be twice as fast, but otherwise unchanged. This is confirmed in figure [7J 
where entanglement transport for various transition dipole strengths is shown, each 
scaled by (|34|). 

Consider next the distance dependence. Numerical calculations are shown in 
figure [HI Here we have kept the ratio a/xo constant and scaled the distances between 
the atoms according to (I3ip . however we did not scale the the width o"o of the initial 
nuclear wavefunction. One sees that the overall dynamics obeys the scaling ( 134]) . however 
there are slight differences in the magnitude of the entanglement. These are due to the 
different relative width of the initial Gaussian. If we scale also <7o we obtain perfect 
agreement as in figure [7J 

The dependence on the mass is demonstrated in figure [HI Again one sees that 
the scaling relation ([34]) is basically fulfilled, however, there are quite large changes 
in the absolute value of entanglement when masses or calculation methods are varied, 
especially after the delocalisation of the excitation (see section |3~2|) . These differences 
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Figure 8. As figure[7]but now for scaled distances A x 2/mi and constant cl/xq. Here 
the width oq of the initial nuclear wave- function is kept constant. The time is given 
in units of t = t (X) = TA 5 / 2 with T = 6.44^s. 



appear pronounced for very small massegj. 

To understand the deviations, the number of jumps invoked in Tully's method is 
shown as a function of mass in figure [101 As expected for smaller masses the dynamics 
becomes more non- adiabatic and the number of jumps strongly increase iffl . 

From the above analysis one sees that due to an increasing number of jumps 
there are differences in the dynamics for different masses. However, in practice these 
differences are quite small as is demonstrated in figure [9j The overall shape of the curves 
stays the same, only the timescales for the total dynamics changes. Nevertheless one 
clearly sees that with smaller mass, where the dynamics becomes more non-adiabatic, 
there is also a larger loss of entanglement. 

It is instructive also to compare Tully's surface hopping approach with the much 
simpler (and faster) Ehrenfest method. For larger masses, in figure [9j the results 
obtained from Tully's surface hopping method and the Ehrenfest-average-potential 
method do not differ, since only a few transitions between the adiabatic states occur. 
For lighter masses, where more transitions occur, the methods start to deviate, since the 

jj The green curve corresponds roughly to hydrogen. 

ff A similar plot but with the number of jumps reduced by roughly a factor of ten is obtained 
when starting from the symmetric initial electronic state, due to a larger energy-separation between 
neighbouring surfaces. 
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Figure 9. As figure[7]but now for different masses M = /3Mu. Tully's surface hopping 
algorithm (solid lines), Ehrenfest method (dashed lines). The time is given in units 
of the mass-dependent time to = to((3) = TyfP with T = 6.44yus. For m = 1 and 
m = 1.49 all lines are indistiguishable. 



dynamics on the averaged potential in the Ehrenfest-method differs from that stemming 
from Tully's algorithm. 

4-1. Dependence on the interaction potential 

The physical transition dipole-dipole interaction between the states | ir n ) scales with 
1/R 3 , where R is the distance between two Rydberg atoms. In the following we 
will investigate (hypothetical) resonant energy transfer interactions with power law 
dependence on the distance 

V nm = -ii 2 /R a nm . (35) 

As already shown in figure [2l to obtain the initial state | ipi n i ) ~ (| ix\ ) — \ t\2 ))/ \/2, the 
ratio a/x has to be decreased for decreasing a. To fix the a dependence of a/x by 
keeping the ratio of the interaction energy between atom 1 and 2 to that between atom 
2 and 3 the constant compared to the case a = 3. Parameters selected in this manner 
are shown in figure [2] as solid black bars. 

To have comparable dynamics for different exponents a, we further fix the values of 
a and x . We choose the a dependent initial distance a a by requiring Vi^aQ.) = ^12(03), 
which leads to a a = (a-s) 3 ^ a . Furthermore we have to adapt the width of the initial 
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Figure 10. Number of jumps invoked in Tully's algorithm as function of the mass. 
Here j3 — M/Mn is the mass expressed in units of the mass of Lithium. Black squares 
are the numerical data (with statistical errors) and the red lines are power law fits. 



nuclear wave packet to the new distance. Our variance a a of the Gaussian position 
distribution is determined by requiring that the corresponding spread of initial potential 
energies is roughly independent of a. This leads to the condition a a ~ - a% cr 3 . 

In figure [11] the entanglement between neighbouring atoms is shown for different a 
and M = Mia- The time is in units of to(a), which is different for each a 

to(a) = - y/fadaY-Z T (36) 

with with T = 6.44/is. One sees that in these scaled time units, for a > 3, the 
transport of entanglement is more or less independent of a. For smaller a, however, the 
entanglement transport is strongly reduced. 

5. Conclusions and Outlook 

We have demonstrated a strong connection between the motion of a chain (aggregate) 
of Rydberg atoms and the coherent propagation of a single electronic excitation within 
the chain. Adiabatic transport ensures that the excitation remains spatially localised 
near a dislocation passing through the chain. Our results were obtained with Tully's 
surface hopping method [5U [55] and the Ehrenfest method [531 Ej, both of which we 
vindicate by comparison with exact quantum calculations for a smaller model system 
with similar dynamics. A key feature of our setup is that the initial state is a repulsive 
electronic eigenstate of the chain. 

If the system was not prepared in an electronic eigenstate but in a state where the 
excitation is localised on a single atom, one would find a fast excitation transfer similar 
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Figure 11. As figure[7]but now for different exponents a. The time is scaled according 

to dSU). 



to that described e.g. in [2U [58]. If free motion is added to such a scenario, however, 
those parts of the population that necessarily initially reside on an an attractive surface 
can lead to fast collisions of light atoms. In that case the dynamics could be treated 
with the Tully algorithm, however the Ehrenfest method would fail (e.g. for the dimer 
there would be no movement at all). 

We wish to contrast the results of the present study with those obtained in [10] 
for van-der-Waals interactions. For the case of repulsive van-der-Waals interaction one 
finds similar trajectories for the movement of the atoms, showing for example Newton's 
cradle like transfer of a dislocation through an atomic chain. The crucial difference to 
the dipole-dipole coupling presented here is the excitation energy transfer involved in 
the latter. Even more important is that the dynamics in the dipole-dipole case depends 
on the electronic state. An even stronger contrast to van-der-Waals is found for initial 
states other than the fully repulsive one treated here. Even mixed, partially attractive, 
partially repulsive dynamics is possible as shown by Ates et al. [TS]. Such dynamics 
would arise from scenarios presented here only if a sufficient fraction of the atoms has 
undergone non-adiabatic transitions to other potential surfaces (as in our figure HJd)). 

In this article we have exclusively studied free Rydberg atoms. When the atoms are 
trapped, the dipole-dipole forces will induce oscillations of atoms in the traps, which 
in turn lead to oscillating couplings and again to a correlation between the motion 
of the exciton and the motion of the Rydberg atoms. This gives rise to the well 
known Davydov-soliton [59l [60]. Constructing large "crystals" of N Rydberg atoms 
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is problematic since each atom has a finite life-time r, leading to an even shorter life- 
time of the crystal of r/N. For the parameters used in our simulations of section [31 
we expect the lifetime to be sufficient. Therefore it would be advantageous to map 
the strong dipole-dipole interaction in the Rydberg state to the ground state using 
off-resonant laser dressing techniques [61]. Using this technique and a ring-geometry, 
it is even possible to use (dressed) Rydberg aggregates for the study of near conical- 
intersection dynamics [62J. 

Experimentally, observables as shown in figure [6] could be monitored using 
techniques for the simultaneous position and state measurement of Rydberg atoms [63J . 
In our particular system the presence of entanglement can then be directly inferred from 
the state populations. 



Appendix A. Entanglement of formation 

The "entanglement of formation" [571 [64] is an entanglement measure for bipartite 
quantum states, also applicable to mixed states. For a pure state it equals 1 for perfect 
entanglement and it is for a separable state. 

We calculate this entanglement measure in the following way: First consider the 
reduced density matrix describing the electronic state when the atomic positions are 
traced out 



a 



/s&nml K n ) { lT m | , (A.l) 



with 



/ 0*(R)0 m (R) QM 
o- nm = { (A.2) 
c* c m Tully / Ehrenf est . 

The first expression holds for the full quantum calculations, the second for the quantum- 
classical methods. In the latter case 7TT denotes the trajectory average and c n = 
Ylk^nk^k are the coefficients in the diabatic basis with O nfc defined in ( I2"2"j) . 

From flA.2j) we construct the binary reduced electronic density matrix of atoms a 



and b 

P ab = Tr^[a]. (A.3) 

The symbol Tr^ a '^ [ • • • ] denotes the trace over the electronic states for all atoms other 
than a, b. Recall that in the present approach each atom is described by a two level 
system as discussed in section 12.21 to 12.61 With our labels for those two states, | s ) 
and | p ) , the trace appearing in ( 1A.3|) is over the Hilbert space spanned by the basis 



{\n 1 )---\n N ), rij G s,p} 

The remaining reduced subspace of atoms a and b is spanned by \pp), \ps), \ sp), 
\ss). Since the states | n n ) appearing in ( lA.ljl only contain a single excitation p, all 
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matrix elements of f3 ab involving | pp ) vanish. In the reduced basis one finally has 

/ \ 

(J aa a ab 
a* ab a bb 

\ E^w> } ^7 

From this we construct 



oh 



(A.4) 



S a b — \l \/ PabPabV flab, 

C ab = max(0, 2X ab - TrS ab ), 
With the further definitions 

H(x) — — [x log 2 x + (1 — x) log 2 (l — x)], 



(A.5) 
(A.6) 

(A.7) 
(A.8) 



£{x) = H(l/2 + Vl-x 2 /2), 

where X ab denotes the largest eigenvalue of S ab , we can finally obtain the binary 
entanglement of the electronic states of atoms a and b 

E ab = S(C ab ). (A.9) 

For further details we refer to [51 



Appendix B. Tully's surface hopping 

The quantum mechanical dynamics governed by (124")) and the classical equation of 
motion (1271) are solved self-consistently. The atoms move on a single adabatic potential 
surface k, which however may be changed via sudden jumps to another surface q. The 
probability for a jump from state k to state q is given by 



g kq = max(0 



(B.1) 



where At denotes the propagation time step and 

b ak =-2Re(a* qk R-d qk ), (B.2) 

a qk = C q C* k . (B.3) 

To determine if during a time step a jump takes place we compare g kq with a uniform 
random number £ £ [0, 1]. If £ < g kl , the jump is to the surface q — 1, if 
9ki < £ < gki + gk2, to q = 2 and so forth. When a switch takes place the velocity 
R is adjusted in order to conserve the total amount of energy. This will be done in the 
direction of the non-adiabatic coupling vector d kq as follows 

Ikqd-kq 



Kit) = R(t - At) 



M 



Here 



(3 kq ±J(3l Q + Aa kq (U k -U q 



Ikq 



2a 



kq 



Pkq^O 



(B.4) 



(B.5) 
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with 

1 N 

AT 

fl^X^-dg. (B.7) 

1=1 

If the energy of the final surface q is larger than that of the initial surface k and the 
velocity reduction required is greater than the component of velocity to be adjusted, 
then the jump is rejected and instead of (IB.5P we use 

l kq = ^- (B.8) 

Ukq 

which corresponds to a reflection of the velocity component along d^. Further details 
about Tally's algorithm can be found e.g. in [5U 155] . 

As described in [65] the forces V R C/ fe (R) = V R (^ fc (R)|# cl (R)|^ fc (R)) and the non- 
adiabatic coupling vectors are calculated using the Hellman-Feynman theorem. One 
then finds 

V R f/ fe (R) = (^ m (R)|(V R i/ cl (R))|^ m (R)). (B.9) 

and 

(^(R)|(V R ff cl (R))|^-(R)) 

dfc9 = — urn-urn — (R10) 

The total density n(R, t), as defined in (|23p . is obtained through a binning of the 
single trajectories, which means for each time step our spatial domain is discretised into 
bins and if a trajectory R(t) = (R!(t)...R N (t)) T lies within such a bin, n(R,t) for that 
bin will be increased by one. By normalising n(R, t) one obtains the probability to find 
a atom at a given time in a certain interval of space. 



Appendix C. Phenomenological model of ionisation 

Our essential states model as justified in section 12.41 is only valid while atoms do not 
approach each other closely. Once they do, dipole-dipole shifts of all electronic states 
become too large to work in a small Hilbert-space of electronic states. The most 
prominent consequence is collisional ionisation of Rydberg atoms [10]. In order to avoid 
excursions of our numerical propagation schemes into realms where the underlying model 
is invalid, we incorporate a simple, phenomenological treatment of ionisation for very 
close atoms. 

For the quantum mechanical calculations, we employ an imaginary absorbing 
potential of the form 

i J^n(R) = • • • - iW(R)0„(R). (C.l) 

into (|20|) . The shape of VT(R) is chosen to minimise reflection, while fully removing 
components of the wave function that correspond to atoms closer than an enforced 
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minimal distance. For the data of figure Hp-j), where the ionisation is most important, 
this distance is 2/xm. For the two quantum-classical trajectory methods (EF, Tally), we 
incorporate the effect described by flClH through a stochastic "ionisation probability" 
2W(x)At in each discrete time-step of duration At. 

Note that we do not aim to model realistic ionisation rates, however we do employ 
the same model of ionisation in all three methods (QM, EF, Tally). Further we point out 
that this approach practically leads to ionisation of all iV atoms, even if only two atoms 
collided. A physically correct treatment would require a density matrix formalism, going 
far beyond our goal, to simply to exclude numerical data from unrealistic regions of the 
model. 



References 

[1] M. D. Lukin, M. Fleischhauer, R. Cote, L. M. Duan, D. Jaksch, J. I. Cirac and P. Zoller; Dipole 
Blockade and Quantum Information Processing in Mesoscopic Atomic Ensembles; Phys. Rev. 
Lett. 87 037901 (2001). 

[2] E. Urban, T. A. Johnson, T. Henage, L. Isenhower, D. D. Yavuz, T. G. Walker and M. Saffmann; 

Observation of Rydberg blockade between two atoms; Nature Physics 5 110 (2009). 
[3] A. Gaetan, Y. Miroshnychenko, T. Wilk, A. Chotia, M. Viteau, D. Comparat, P. Pillet, 

A. Browaeys and P. Grangier; Observation of collective excitation of two individual atoms in the 

Rydberg blockade regime; Nature Physics 5 115 (2009). 
[4] C. H. Greene, A. S. Dickinson and H. R. Sadeghpour; Creation of Polar and Nonpolar Ultra-Long- 
Range Rydberg Molecules; Phys. Rev. Lett. 85 2458 (2000). 
[5] I. C. H. Liu, J. Stanojevic and J. M. Rost; Ultra-Long-Range Rydberg Trimers with a Repulsive 

Two-Body Interaction; Phys. Rev. Lett. 102 173001 (2009). 
[6] V. Bendkowsky, B. Butscher, J. Nipper, J. B. Balcwski, J. P. Shaffer, R. Low, T. Pfau, W. Li, 

J. Stanojevic, T. Pohl and J. M. Rost; Rydberg Trimers and Excited Dimers Bound by Internal 

Quantum Reflection; Phys. Rev. Lett. 105 163201 (2010). 
[7] W. R. Anderson, J. R. Veale and T. F. Gallagher; Resonant Dipole-Dipole Energy Transfer in a 

Nearly Frozen Rydberg Gas; Phys. Rev. Lett. 80 249 (1998). 
[8] W. Li, P. J. Tanner and T. F. Gallagher; Dipole-Dipole Excitation and Ionization in an Ultracold 

Gas of Rydberg Atoms; Phys. Rev. Lett. 94 173001 (2005). 
[9] F. Robicheaux, J. V. Hernandez, T. Topcu and L. D. Noordam; Simulation of coherent interactions 

between Rydberg atoms; Phys. Rev. A 70 042703 (2004). 
[10] T. Amthor, M. Reetz-Lamour, C. Giese and M. Weidcmiiller; Modeling many-particle mechanical 

effects of an interacting Rydberg gas; Phys. Rev. A 76 054702 (2007). 
[11] H. Schcmpp, G. Giintcr, C. S. Hofmann, C. Giese, S. D. Saliba, B. D. DePaola, T. Amthor, 

M. Weidemuller, S. Scvincli and T. Pohl; Coherent Population Trapping with Controlled 

Interparticle Interactions; Phys. Rev. Lett. 104 173602 (2010). 
[12] I. Lesanovsky; Many-Body Spin Interactions and the Ground State of a Dense Rydberg Lattice 

Gas; Phys. Rev. Lett. 106 025301 (2011). 
[13] N. Henkel, R. Nath and T. Pohl; Three- Dimensional Roton Excitations and Supersolid Formation 

in Rydberg- Excited Bose-Einstein Condensates; Phys. Rev. Lett. 104 053004 (2010). 
[14] O. Miilkcn, A. Blumen, T. Amthor, C. Giese, M. Reetz-Lamour and M. Weidcmiiller; Survival 

probabilities in coherent exciton transfer with trapping; Phys. Rev. Lett. 99 090601 (2007). 
[15] C. Ates, A. Eisfeld and J. M. Rost; Motion of Rydberg atoms induced by resonant dipole-dipole 

interactions; New J. Phys. 10 045030 (2008). 



Adiabatic entanglement transport in Rydberg aggregates 28 

[16] J. Frenkel; Zur Theorie der Resonanzverbreiterung von Spektrallinien; Zeitschrift fur Physik A 59 
198 (1930). 

[17] J. Frenkel; On the Transformation of Light into Heat in Solids. I; Phys. Rev. 37 17 (1931). 

[18] J. Franck and E. Teller; Migration and Photochemical Action of Excitation Energy in Crystals; J. 

Chcm. Phys. 6 861 (1938). 
[19] A. Davydov; Theory of Molecular Excitons; McGraw-Hill (1962). 

[20] H. van Amerongcn, L. Valkunas and R. van Grondcllc; Photo synthetic Excitons; World Scientific, 

Singapore (2000). 
[21] T. Kobayashi, editor; J- Aggregates; World Scientific (1996). 

[22] A. Bierman; Motion of Vibrationally Coupled Excition Wave Packet in an Infinite 1-Dimensional 

Crystal; J Chcm Phys 46 1484 (1967). 
[23] H. Haken and P. Reineker; The coupled coherent and incoherent motion of excitons and its influence 

on the line shape of optical absorption; Z. Phys. 249 253 (1971). 
[24] J. Roden, A. Eisfeld, W. Wolff and W. T. Strunz; Influence of Complex Exciton-Phonon Coupling 

on Optical Absorption and Energy Transfer of Quantum Aggregates; Phys. Rev. Lett. 103 058301 

(2009). 

[25] P. Rebentrost, R. Chakraborty and A. Aspuru-Guzik; Non-Markovian quantum jumps in excitonic 
energy transfer; J. Chem. Phys. 131 184102 (2009). 

[26] A. Blumen and R. Silbey; Exciton line shapes and migration with stochastic exciton lattice coupling; 
J. Chem. Phys. 69 3589 (1978). 

[27] I. Barvi'k and V. Szocs; Coherent and incoherent exciton motion in the framework of the continuous- 
time random walk; Physics Letters A 125 339 (1987). 

[28] F. Hofelich; Die Bewegung eines Exzitons entlang eines Polymers unter dem Einflu der 
Gitterschwingungen; Z. Phys. B: Physik Der Kondensierten Materie 5 208 (1966). 

[29] A. Eisfeld and J. S. Briggs; The J-band of organic dyes: lineshape and coherence length; Chem. 
Phys. 281 61 (2002). 

[30] F. Caruso, A. W. Chin, A. Datta, S. F. Huelga and M. B. Plcnio; Highly efficient energy excitation 

transfer in light-harvesting complexes: The fundamental role of noise-assisted transport; The 

Journal of Chemical Physics 131 105106 (2009). 
[31] M. Sarovar, A. Ishizaki, G. R. Fleming and K. B. Whaley; Quantum entanglement in photo synthetic 

light harvesting complexes; Nature Physics 6 462 (2010). 
[32] G. S. Engel, T. R. Calhoun, E. L. Read, T.-K. Ahn, T. Mancal, Y.-C. Cheng, R. E. Blankenship 

and G. R. Fleming; Evidence for wavelike energy transfer through quantum coherence in 

photo synthetic systems; Nature 446 782 (2007). 
[33] H. Lee, Y.-C. Cheng and G. R. Fleming; Coherence dynamics in photosynthesis: protein protection 

of excitonic coherence; Science 316 1462 (2007). 
[34] E. Collini and G. D. Scholcs; Coherent intrachain energy migration in a conjugated polymer at 

room temperature; Science 323 369 (2009). 
[35] W. H. Li, P. J. Tanner and T. F. Gallagher; Dipole-dipole excitation and ionization in an ultracold 

gas of Rydberg atoms; Phys. Rev. Lett. 94 173001 (2005). 
[36] S. Wiister, C. Ates, A. Eisfeld and J. M. Rost; Newton's cradle and entanglement transport in 

flexible Rydberg aggregates; Phys. Rev. Lett. 105 053004 (2010). 
[37] A. Asadian, M. Tiersch, G. G. Guerreschi, J. Cai, S. Popescu and H. J. Bricgel; Motional effects 

on the efficiency of excitation transfer; New J. Phys. 12 075019 (2010). 
[38] M. Schwocrcr and H. Wolf; Organic Molecular Solids; Wiley- VCH (2006). 

[39] D. M. Eisele, J. Knocster, S. Kirstein, J. P. Rabe and D. A. Vanden Bout; Uniform exciton 
fluorescence from individual molecular nanotubes immobilized on solid substrates; Nature 
Nanotcchnology 4 658 (2009). 

[40] T. Tani; J- Aggregates in Spectral Sensitization of Photographic Materials; in T. Kobayashi, editor, 
J-Aggregates; World Scientific (1996). 

[41] K. Takechi, P. K. Sudcep and P. V. Kamat; Harvesting infrared photons with tricarbocyanine dye 



Adiabatic entanglement transport in Rydberg aggregates 



29 



clusters; J Phys Chem B 110 16169 (2006). 
M. Recrs, T. W. Smith and L. B. Chen; J-aggregate formation of a carbocyanine as a quantitative 

fluorescent indicator of membrane-potential; Biochemistry 30 4480 (1991). 
V. N. Dedov, G. C. Cox and B. D. Roufogalis; Visualisation of mitochondria in living neurons 

with single- and two-photon fluorescence laser microscopy; Micron 32 653 (2001). 
S. N. Rylova, A. Amalfitano, D.-A. Persaud-Sawin, W.-X. Guo, J. Chang, P. J. Jansen, A. D. Proia 

and R.-M. Boustany; The CLN3 Gene is a Novel Molecular Target for Cancer Drug Discovery; 

Cancer Research 62 801 (2002). 
S. Dahne; Nanostrukturierte J-Aggregate als potenzielle Lichtsammelsysteme fur Photosynthesen; 

Bunsen-Magazin 4 81 (2002). 
S. Kirstein and S. Dahne; J-aggregates of amphiphilic cyanine dyes: Self- organization of artificial 

light harvesting complexes; International Journal of Photoenergy page 20363 (2006). 
J. Rodcn, A. Eisfeld, M. Dvorak, O. Biincrmann and F. Stienkemeier; Vibronic line shapes of 

PTCDA oligomers in helium nanodroplets; Journal of Chemical Physics 134 054907 (2011). 
G. S. Engel, T. R. Calhoun, E. L. Read, T. K. Ahn, T. Mancal, Y. C. Cheng, R. E. Blankenship 

and G. R. Fleming; Evidence for wavelike energy transfer through quantum coherence in 

photosynthetic systems; Nature 446 782 (2007). 
G. Birkl and J. Fortagh; Micro traps for quantum informa- tion processing and presicion force 

sensing; Laser & Photon. Rev. 1 12 (2007). 
G. D. Stevens, C. Iu, T. Bergeman, H. J. Metcalf, I. Seipp, K. T. Taylor and D. Dclande; Precision 

measurements on lithium atoms in an electric field compared with R-matrix and other Stark 

theories; Phys. Rev. A 53 1349 (1996). 
T. Amthor, M. Reetz-Lamour, S. Westermann, J. Denskat and M. Weidemuller; Mechanical Effect 

of van der Waals Interactions Observed in Real Time in an Ultracold Rydberg Gas; Physical 

Review Letters 98 023004 (2007). 
J. Kruse, C. Gierl, M. Schlosser and G. Birkl; Reconfigurable site-selective manipulation of atomic 

quantum systems in two-dimensional arrays of dipole traps; Phys. Rev. A 81 060308(R) (2010). 
M. S. Topaler, T. C. Allison, D. W. Schwenke and D. G. Truhlar; What is the best semiclassical 

method for photochemical dynamics of systems with conical intersections?; J. Chem. Phys. 109 

3321 (1998). 

J. C. Tully and R. K. Preston; Trajectory Surface Hopping Approach to Nonadiabatic Molecular 
Collisions: The Reaction of H+ with D 2 ; J. Chem. Phys. 55 562 (1971). 

J. C. Tully; Molecular dynamics with electronic transitions; J. Chem. Phys. 93 1061 (1990). 

M. A. Nielson and I. L. Chuang; Quantum computation and quantum information; Cambridge 
University Press (2000). 

S. Hill and W. K. Wootters; Entanglement of a Pair of Quantum Bits; Phys. Rev. Lett. 78 5022 
(1997). 

A. Eisfeld; Phase directed excitonic transport and its limitations due to environmental influence; 

Chemical Physics 379 33 (2011). 
A. S. Davydov and N. I. Kislukha; Solitary Excitons in One- Dimensional Molecular Chains; phys. 

stat. sol. (b) 59 465 (1973). 
W. Weidlich and W. Hcudorfer; Coherent Exciton-Phonon Pulse Propagation in One Dimension; 

Z. Physik 268 133 (1974). 
S. Wiistcr, C. Atcs, A. Eisfeld and J. M. Rost; Excitation transport through Rydberg dressing; 

|http : //arxiv . org/abs/1011 . 5483| (2010). 
S. Wiistcr, A. Eisfeld and J. M. Rost; Conical intersections in an ultracold gas; Phys. Rev. Lett., 

in press; |http : //arxiv . org/abs/1011 . 5489| (2011). 



C. van Ditzhuijzen, A. Koendcrink, L. Noordam and H. van Linden van den Heuvell; Simultaneous 
position and state measurement of Rydberg atoms; Eur. J. Phys. D 40 13 (2006). 

W. K. Wootters; Entanglement of Formation of an Arbitrary State of Two Qubits; Phys. Rev. 
Lett. 80 2245 (1998). 



Adiabatic entanglement transport in Rydberg aggregates 



30 



[65] S. Hammes-Schiffer and J. C. Tully; Proton- Transfer In Solution - Molecular- Dynamics With 
Quantum Transitions; J. Chem. Phys. 101 4657 (1994). 



